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Abstract. We construct center-stable and center-unstable manifolds, as well 
as stable and unstable manifolds, for the nonlinear Klein-Gordon equation with 
a focusing energy sub-critical nonlinearity, associated with a family of solitary 
waves which is generated from any radial stationary solution by the action of 
all Lorentz transforms and spatial translations. The construction is based on the 
graph transform (or Hadamard) approach, which requires less spectral information 
on the linearized operator, and less decay of the nonlinearity, than the Lyapunov- 
Perron method employed previously in this context. The only assumption on the 
stationary solution is that the kernel of the linearized operator is spanned by its 
spatial derivatives, which is known to hold for the ground states. The main novelty 
of this paper lies with the fact that the graph transform method is carried out in 
the presence of modulation parameters corresponding to the symmetries. 
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1. Introduction 
Consider the focusing nonlinear Klein-Gordon equation (NLKG) on Mf' 

u-Au + u = f{u), u{t,x) ■.R^^'^ (1.1) 

where / : M — )• M is a given nonlinearity. A typical example is the focusing power 
nonlinearity 

f(u) = \u\P-\ 3<p+l<l^ id>3) 

^ ' ' \oo {d<2). ^ ' 

The lower bound can in principle be reduced to p > 1, but we assume p > 2 to avoid 
technical and non-essential complications in the nonlinear estimates. 
The equation preserves the total energy and momentum 

E{u):= [ [y^±r^!^^±y!-/(-i)(n)]da;, P{u) := f iiVudx, (1.3) 

where /^^""^^ : M — )■ M is the primitive /^^""^^(a) = f{b)db. These quantities are 
well-defined in the energy space 

u{t) := {u{t),u{t)) eH:= H\W^) x L^R"^). (1.4) 

Throughout the paper, we do not distinguish vertical and horizontal vectors in T-L, 
unless it may lead to any confusion. 

We will consider NLKG in the energy space Ti, regarding it as a Hamiltonian 
system. Our goal is to construct a local center-stable manifold of the family of 
traveling waves generated by the Lorentz transforms and spatial translations acting 
on a stationary solution. For brevity, we call the latter manifold of traveling waves 
the soliton manifold. There are two major approaches used in the construction of 
center-stable manifolds: the Hadamard method and the Lyapunov-Perron method. 
The former uses the evolution backward and locally in time to find a flow-invariant 
graph of the unstable modes in terms of the other components (the Hadamard 
approach also goes by the name of graph transform or invariant cones method). The 
latter uses the evolution forward and globally in time to find an initial adjustment 
of the unstable modes so that they remain small forever. 

Bates and Jones [Ij developed the Hadamard method in the general setting of 
an ODE of the form x = Ax + f{x) where A is an (unbounded) operator on some 
Banach space X which generates a continuous semi-group, the nonlinearity / is 
locally Lipschitz on X, satisfies /(O) = 0, and admits arbitrarily small Lipschitz 
constants near the equilibrium x = 0. The spectrum of A is divided in the stable 
part with eigenvalues in the left-half plane, the unstable part which lies in the right- 
half plane, and the center part which lies on the imaginary axis. Assumptions are 
made on the dimensions of the corresponding spectral subspaces of X, and the 
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associated flows (if tfie spaces are infinite-dimensional) so as to represent two main 
scenarios: the dissipative case (D) on the one hand, and the conservative case (C) 
on the other hand. For (D) one demands that only the stable subspace be infinite- 
dimensional and that the associated semigroup is exponentially stable. For (C) 
only the center subspace is infinite dimensional, which is precisely what occurs in 
Hamiltonian problems. 

Bates and Jones then verified that the abstract center-stable manifold which they 
constructed in [1] applies to stationary solutions of NLKG under the radial symmetry 
restriction, for the power nonlinearity (11. 2p with p < d > 3, where the upper 
bound on p was required to ensure that the nonlinearity is locally Lipschitz — > L^. 
They also showed that if the linearized operator has no nonzero radial functions in 
its kernel, then 

(1) Every solution starting on the center-stable manifold stays there forever in 
positive times, remaining in a small neighborhood of the stationary solution. 

(2) Every solution starting in that small neighborhood, but off the manifold, 
must exit the neighborhood in finite positive time. 

The kernel condition holds for the ground state (the positive stationary solution), 
by work of Weinstein [H] . 

Gesztesy, Jones, Latushkin, and Stanislavova [8] demonstrated that the Bates- 
Jones theory applies to the nonlinear Schrodinger equation (NLS) with a spatially 
localized nonlinearity. Notice that the radial restriction for NLKG prohibits both 
the spatial translations and the Lorentz transforms, and so the soliton manifold is 
reduced to a fixed stationary solution. Similarly, the localized nonlinearity of [8] 
destroys the scaling, translation, and Galilean invariance, so that the soliton can 
change only with respect to the phase parameter. Indeed, as we will explain below, 
moving sohtons represent a serious obstacle to the Bates- Jones approach. 

On the other hand, the second author [13] developed the Lyapunov-Perron (LP) 
method for the ground state of the cubic NLS in M^, without imposing any symmetry 
assumptions, but in a weighted if**-space (or an unweighted L^-based space). In this 
approach, the soliton is allowed to move. Recently, Beceanu [3] extended the latter 
work to the critical space H^^^ which is bigger than the energy space. Finally, 
and partly based on a novel approach to linear dispersive estimates developed by 
Beceanu P], the authors proved in [12] that the LP approach can be carried out for 
NLKG in the energy space without any symmetry restrictions. 

However, an essential difficulty in applying the LP method to a nonlinear disper- 
sive equation (without dissipation) is that it requires global dispersive estimates, 
which in turn necessitates fine spectral information, such as the absence of threshold 
resonances and of so-called spuriou^ eigenvalues; alternatively, in the presence of 
such spurious eigenvalues one might hope to invoke the Fermi golden rule. Those 
conditions are in general very hard to check, even for the ground state (apart from 
the one- dimensional case [9j where purely analytical arguments are available), al- 
though there has been some recent progress in this direction [HI EJ [101 [S]- 
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While the LP method requires stronger ingredients, it also leads to more detailed 
conclusions. More specifically, one obtains that solutions starting on the center- 
stable manifold scatter to the soliton manifold in forward time. In other words, 
the distinction between the LP approach and the Hadamard approach is roughly 
tantamount to the distinction between asymptotic and orbital stability of solitary 
waves. 

In this paper, we employ the Hadamard method in the nonradial setting. Our 
main challenge is to extend the result by Bates and Jones to the family of traveling 
waves, rather than stationary ones. We therefore have to investigate the dynamics 
along the soliton manifold as well, which is usually called the modulational analysis 
in the stability problem of solitons. In our setting, the soliton manifold has 2d 
dimensions corresponding to the relativistic momentum and position vectors in W^. 

Those parameters can be fixed by means of a Lorentz transform which reduces 
the total momentum to zero, and by using a coordinate moving with the soliton. 
In doing so, we encounter a derivative loss due to the translation, i.e., a transport 
term, in the modulated equation for the difference of two solutions, which disables 
the contraction argument for the graphs in the energy norm. This difficulty is not an 
artifact of the coordinate choice, but a natural consequence of the two facts that the 
solitons are translated by the fiow, while the translation is not Lipschitz continuous 
in any Sobolev space. The same problem occurs for any other continuous group 
action involving a coordinate change, such as scaling or rotation. 

We overcome this difficulty by introducing a nonlinear quasi- distance in the energy 
space, for which the spatial translation becomes Lipschitz continuous, while the 
topology remains the same. Using the contraction mapping principle with this 
distance for the continuous spectral part, we are able to carry out the Hadamard 
method in the presence of modulational parameters. 

A more technical issue concerns allowing nonlinearities all the way up to the 
critical power, i.e., for p < {d + 2)/{d — 2), while Bates and Jones assumed 
p < d/{d — 2). This is easily resolved by using the Strichartz estimate for the free 
Klein-Gordon equation, and by relaxing some fiow-invariance conditions by constant 
multiples. 

Since the description of dynamics around the manifolds ((l)-(2) above) is also 
extended to the current setting, we can easily observe that the maximal backward 
evolution of the center-stable manifold is identical, in a small neighborhood, to the 
forward trapping set 7+: the collection of initial data for which the solution (of the 
original NLKG) stays in the small neighborhood for large times. 

In the special case where the soliton manifold is generated from the ground state, 
we can combine the above result with the one-pass theorem in [12] as well as the 
openness and connectedness of the forward scattering set 5+ and the forward blow- 
up set B+, thereby concluding that 7+ separates locally and globally all the solutions 
with energy at most slightly above the ground state energy into iS+ and There- 
fore, the conclusion of [12] is extended to the range 



deN, l + -<p<l + j-^, p>2, (1.5) 
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except for the following scattering statement on 7+: all solutions in 7+ scatter to 
the soliton manifold as t — ?■ oo. This statement was proved in [12] for d = p = 3 hj 
means of the LP method using the following gap property of the linearized operator 
L+- 

(0, 1) n cr(L+) = 0, and there is no threshold eigenvalue or resonance 

This is proved (at least for the radial case) in [5j. Note that the numerical analysis 
of [B] suggests that the absence of threshold resonances and spurious eigenvalues 
fails for some powers in (1 + 4/3, 3) for d = 3, where the LP method without any 
hypothesis (typically the Fermi golden rule) is not so far available. Also note that 
the lower bound l + A/dis required by the proof of the one-pass theorem, but not by 
the Hadamard construction in this paper, while the Lyapunov-Perron method also 
needs it in order to work in the energy space. 

To state the main result, we clarify the assumptions on the nonlinearity / and on 
the stationary solution: 

feC\R;R), = /(O) = /'(O), 

(l + \ar' (rf>3,2<3p<^) 
VaGM, \f"{a)\< ll + \a\P~^ (rf = 2, 2 < 3p < oo) 
I arbitrary (d = 1). 

To have manifolds, we assume a bit more regularity: for some p > 2 in the above 
range, 

ir(ai) - r(a2)| < (|ai - aal^-' + ki - a2|)[l + {a^r' + \a2r\ (1.7) 
These assumptions are satisfied for example by 

A;;finite 

provided that all pk > 2 are in the range fll.6p . 

Let Q G H\R'^) be a stationary solution of NLKG, i.e., a weak solution of the 
elliptic PDE 

-AQ + Q = f{Q). (1.9) 

Standard arguments imply that Q G with exponential decay as |x| — )■ oo. 
For existence, see the classical work by Berestycki, Lions P]. The action of the 
Lorentz transforms and the spatial translations generate the traveling wave family 
parametrized by the relativistic momentum p G M'' and position q G M*^: 

Q{p,q) ■=Q{x-q + p{{p)-l)\p\-^p-{x-q)), (1.10) 

so that each traveling wave can be written in the form 

u{t) = Q{p,q{t)). Jt^(t) = ^y (1-11) 
For brevity, the spatial translate is denoted also as 

Qcix) ■.= Q{x-c). (1.12) 
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The vector form is denoted by 

Q := (Q, 0), Q(p, q) ■= (Q(p, ■ Vg(p, g)), (1.13) 

and the sohton manifold of Q is defined as 

y{Q) ■■= {QiP, ?1 WeM^ c n, (1.14) 

which is a manifold of dimension 2d. The linearized operator at Q 

L^:=D'-f'iQ) = -A + l-fiQ), D := (1.15) 

is self-adjoint on with a finite number of eigenvalues and continuous spectrum 
ac{L^) = (Tac{L^) = [l,oo). The translation invariance of NLKG implies that 
Lj^VQ = 0. The only assumption on Q in this paper is 

L;1(0) = span{Vg}. (1.16) 

This is a well-known property of the ground states. To be more precise, by an 
argument of Weinstein [14j, it holds for the ground state Q provided no radial 
function lies in the kernel of L+. The latter holds for any subcritical monomial 
nonlinearity (as well as others), see Lemma 2.3 in [11], for example. Moreover, 
( I1.16P seems to be a natural assumption for any other radial static solution. For 
non-radial static solutions, we have to include angular derivatives as well, but we 
do not consider such solutions in this paper. Although we will not explicitly use the 
radial symmetry of Q, the reader may assume it without losing anything throughout 
the paper. 

Theorem 1.1. Let d & N and assume that f satisfies f ll.6p . Let Q be a static 
solution (11. 9p and assume that its linearized operator satisfies fll.l6p . Then 
there is a Lipschitz manifold A4cs containing the soliton manifold S^{Q), with 

the following properties: 

(1) The codimension of M.cs in equals the total dimension of the eigenspaces 
of corresponding to negative eigenvalues, which we denote by K. 

(2) M.CS is invariant under the forward evolution of NLKG (11. ip . 

(3) M.CS is invariant under spatial translations. 

(4) For every Lorentz transform and every {p,q), there is a small neighborhood 
of Q{p, q) such that the Lorentz transform of any forward global solution 
starting from Aics within this neighborhood remains on Aics for all t > 0. 

(5) Aics is normal at Q to the vector {—kp,p) for any p G solving L+p = 
—k'^p for some k > 0. In other words, 

{u2\p) = ~k{ui - Q\p) + o(||mi - QWm + \\u2\\l-^)- (1.17) 

(6) For any neighborhood O of S^{Q), there is a smaller neighborhood O' , such 
that every solution starting from CflA^cs remains in OnM.cs for allt > 0. 

(7) There is a neighborhood O of S^{Q) such that every solution starting from 
0\A4.cs exits O in finite positive time. 

(8) // in addition f satisfies (ll.7p . then Jvl^s is C^'", where a = min(l,p — 2). 
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The corresponding statement for a center-unstable manifold follows simply by the 
time inversion, so we omit it. However, in the proof we will actually consider the 
center-unstable manifold, for which the forward evolution is used as a contraction 
mapping in the Hadamard method. The center manifold is obtained by intersecting 
the center-stable with the center-unstable manifold. It is of codimension 2K and is 
bi- invariant. 

Properties (6), (7) characterize Aics as the set of solutions which stay close to 
^{Q) for all t > 0. Since the Lyapunov- Perron method looks for such solutions from 
the beginning, it will yield a subset of A4cs, and indeed the same manifold (locally), 
provided that the codimension is the same. An advantage of the Lyapunov-Perron 
method is that it implies the scattering to the soliton manifold for the solutions 
on A4cs (cf. pS])- It will be interesting to see what happens when some spectral 
condition breaks down, e.g., if there is a threshold resonance. 

We also obtain a stable and unstable manifold theorem. Recall the definition of 
K from the previous theorem. 

Theorem 1.2. Under the assumptions of the previous theorem, there exist Lipschitz 
manifolds Aiu an-d Ms of dimensions K + d with the following properties: 

(1) Ms is invariant under the forward evolution of the NLKG f ll.ip . 

(2) Ms is invariant under spatial translations. 

(3) every solution starting on Ms converges exponentially to Q{-—c(t)) as t oo 
where c{t) — )■ exponentially as t oo. 

(4) there exists 6 > small such that Ms is a Lipschitz graph over Bs{0) x M'^, 
with Bs{0) being a 6-ball in the eigenspace corresponding to the negative 
eigenvalues, and the W^- component deriving from spatial translations. 

(5) Mu is obtained from Ms by reversing time. 

2. Preliminaries 

Here we fix some notation. For any two elements f °, in a normed space V, the 
ordered pair and their difference are denoted by 

v"- := {v'',v') eV^, \\v''\\v=\\v'^\\v+\\v'\\v, <v'' := - v'' , (2.1) 

respectively. More generally, for any map M : V x W x ■ ■ ■ and elements G V, 
G W,. . . , the mapped pair is denoted by 

M(v^ w^ . . . ) := {M{v\ w^,...), M{v\w\ ...)). (2.2) 

For any i? G M and 6 > 0, the minimum is denoted by 

(RAS) ■.= mm{R,S). (2.3) 

As usual, a ^ b, a > b, and a c:^ b involve implicit multiplicative constants. 

2.1. Equation and spectrum. The energy space "H = xL^ (Z {L'^Y is endowed 
with the usual inner product 

{if, ij)^ := / [V^i ■ Vipi + ifiitfji + ^2^2] dx (2.4) 
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and the duality coupling 

((^1-0):=/ [(pi{x)i/ji{x) + ip2{x)ilj2{x)]dx, (2.5) 
as well as the symplectic form 

oj{(p,ilj) := (J</?|V') = / [ip2{x)ilJi{x) - (pi{x)ip2{x)] dx, (2.6) 

JRd. 

where J is the skew-symmetric matrix 

Let C be the self-adjoint operator on with domain H'^ x 

Its free version is denoted by 

Then the linearized equation around Q = ((5,0) is 

vt = J>Cv. (2.10) 

The spectrum of J£ is given in terms of that of L+: 

a{JC) = ±^/-a{L+) 

Since /'(Q) is bounded and exponentially decreasing, there are < ^ < 1 < A; and 
a finite set K C [k, k] such that 

a(L+) \ [k^, oo) = {0} U {-A;2 | A; e iT}. (2.11) 

With slight abuse of notation, we let K count the multiplicity of each negative 
eigenvalue —k"^ as well. For each k E K, let pk G 5(M'') be an eigenfunction 
satisfying 

L+pk = -k'^Pk, WPkh = 1- (2.12) 
The (generahzed) eigenfunctions of JjC are 

JCgki = ±kgk±, JCVQ = 0, JCJVQ = -VQ, gk± := f ,\ ^ (2.13) 



which satisfy 

uj{d^Q, JdpQ) = H^M) := {d^QldpQ), uj{gk+,gk-) = 1- (2.14) 
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The corresponding symplectic (spectral) decomposition takes the form 
V = ^ Xk±gk± + ■ VQ + ■ ^VQ + 7) 

Afc± := P±kV ■■= oj{v, ifi-fczp), 

:= p^v := H{Q)-^co{v, JVQ), v := P,v := H{Q)-'u{v, -VQ), 

7 := P^v ■= V - ^ >^±kgk± - /i ■ VQ - u ■ JVQ. 

±,k 

We will use the following projections as well: 

v± := P±v := ^ \k±gk±, Vd := PdV := v - P^v, 

k&K 



(2.15) 



(2.16) 



Vq := PqV := VQ ■ fi + JVQ ■ v>o := P>oV := v - v^, 
:= P^+v := 7 + v+, Vo± := Po±v := Vq + v±, 
and the corresponding subspaces 'H± := P±{'H). Fixing a small number 

< K < fc, (2.17) 
we define the energy norm on to be 

ll^lll — '^'^•l^ + I'^l^ + '^^l^l^ + ^^^1^) - ll^llw' (2.18) 

k£K 

where the final equivalence follows from the orthogonality of 71: 

0=(7i|Vg) = (7i|pfc) Vfc, (2.19) 

together with fll.l6p . since (>C7|7) = (L+71I71) + || 72 Hi- 
Let u E C{P, T-L) be a solution (11. ip . then v = (fi, ^2) := u{x + c) — Q solves 

vt = JCv + c-{VQ + Vv) + N{v), N{v) = {0,N{v)), (2.20) 

where N : ^ H^^ carries the superlinear part 

N{v) := f{Q + v^) - f{Q) - f\Q)v, = o{v^). (2.21) 

We remark that V2 7^ Vi unless c = 0. 

The conserved momentum can be rewritten as 

P{u) := uj{u, Vm)/2 = uj{v, VQ) + uj{v, Vv)/2. (2.22) 

By means of a Lorentz transform, we can reduce the dynamics near the soliton 
manifold ^{Q) to the invariant subspace 

:= {m e -H I P(m) = 0}. (2.23) 

Furthermore, we can restrict v to the subspace 

H^:={veH\ uiy, VQ) + w(t;, Wv)/2 = 0, uj{v, JVQ) = 0}, (2.24) 

by choosing c so that 

= dtco{v, JVQ) = {H{Q) - {V^Q\vi))c + uj{v, VQ). (2.25) 
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Hence the evolution for small v on T-L± is given by 

vt = JCv + A{v) ■ V(g + v) + N{v), 

A{v):={H{Q)-{V'Q\v^))-'co{v,Vv)/2. 

This is a first-order autonomous equation in "H with the superlinear term 

M{v) := A{v)-V{Q + v) + N{v). (2.27) 

In order to implement the Hadamard method, we need to localize the nonlinear 
part M{v) near so that it becomes a small Lipschitz term globally in the energy 
space %. It seems extremely hard to do this keeping the above orthogonal structure, 
since the acceleration or the modulation term is naturally unbounded, unless the 
linearized operator is modified depending on the distance of v from 0. Therefore we 
will not enforce the orthogonality conditions, but instead solve a localized version 
of the above autonomous equation in the whole energy space "H. After constructing 
a center-unstable manifold by the Hadamard method for the localized equation, we 
will restrict that manifold to the subspace "H^ in a small neighborhood of to obtain 
a center-unstable manifold for the true equation. 

In the case of the unstable manifold, the exponential decay of f as t — t- — oo 
ensures that the manifold for the localized equation around falls into 'H±, so that 
we can automatically get the manifold of the true equation. 

2.2. Mobile distance. The most serious obstacle to carrying out the graph trans- 
form method in the nonradial setting results from the contraction step in the con- 
struction of the invariant graphs, where the presence of the unbounded translation 
term causes problems. To remedy this, we introduce the mobile distance on %. 
Heuristically speaking, the standard or Sobolev-type norm is too tight for "hor- 
izontal motion" if i— )■ ip{- + a;o) compared with "vertical motion" ip ^ \ip. The 
mobile distance makes translation just as easy as amplification, without changing 
the topology. 

Definition 2.1. For any continuous increasing function (p : [0, oo) — [0, oo) satisfying 

(1) (j){a) > a. 

(2) a < 26 ^ 0(6) < 4</)(a), 

the mobile distance : Ti x H ^ [0, oo) is defined by 

Mv',v'r := inf \\v'-^ - v\- - q)\\l + klVdl^^lU)'- (2.28) 

Obviously, the infimum in f l2.28p is attained at some q G M'^. m,^ is not really a 
distance, but a quasi- distance on T-L. More precisely, we have 

Proposition 2.2. in f l2.28p is a complete quasi- distance on %, satisfying 

(1) m0(f°,f"'^) > where the equality holds iff = v^. 

(2) ra^{v\v^)=m^{v^,v'). 

(3) m0(f°,f^) < Cd [m<^(f°,f^) + m<^(i;^,f^)] for some absolute constant Cd > 0. 

(4) //m0(f™, f") — 7- (n,m ^ oo) then converges in Ti. 
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Moreover, it satisfies with some absolute constant C > 0, 

\\\v^n - \\v'\\h\ + \\D-\v' - v')\\n < Cm^iv^,v') < CWv' - v'\\n, (2.29) 

where D := a/1 — A. These constants, C and C^, do not depend on the choice of (p. 

Hence defines the same topology as "H, differing only in terms of uniformity. 
For example, for any (p we have 

hm m<^(<^e^"^S -<^e^"^0 = 0(1), (2.30) 
since — ^^e*"^^ = ipe'^^i^i+'^Z^) ^ whereas Hyje*""^^ — (— (/^e*"^!) ||^ = 0{n) unless (p = 0. 

Proof of Proposition [KM (1) and (2) are obvious. For the left-most term of (12.291) . 
and with r^v := f (■ — q), 

IP^Wh - \\v^\\h\ < inf min(||t;° - TqV^Wu, \\v^ - ^g^°||«) 

Q 



(2.31) 



for the second term, 

\\D-\v' - v')\\n < \\D-\v' - r,v')\\n + \\D-\r,v' - v')\\n 

< - Tgv'Wn + \q\\\VD-'v'\\n (2.32) 

<\\v'-Ty\\E + \qm\AE), 

while the right bound in f l2.29p is obvious by choosing q = 0. Next we prove the 
quasi-triangle inequality. For any G "H, there are G M*^ such that 

m^{v'',v^) ~ \\v' - Tq^v^ln + |g^|0(min(||t;°||^, \\v^\n)), (j = 1,2), (2.33) 

since the Ti and E norms are equivalent and the H norm is translation invariant. 
If llf^ll-H ^ ™ii(ll'^^ll'H; il'^^^ilw) then m</,(f °, f-') ~ and so the quasi-triangle 

inequality is obvious. Otherwise, 

M^\v') < \\v' - v\- -q' + q')\\H + \q' - q'\(P{mm{\\v'\\H, \\v'\\n)) 



< _ ^1) _ _ ^2)11^ ^ ^ |g^|0(min(||t;l«, ||t;^||^)) 



(2.34) 



To prove the completeness, let f " be Cauchy in m,^. Then so is D~^v'"' in H by 
(12.29p . Hence, converges to some v G DT-L. (12.291) implies that Hw^H-h converges. 
We may assume that this limit is positive, since otherwise the convergence to is 
obvious. Since f" is bounded in "H, it converges weakly to v in Ti. Passing to a 
subsequence, we may find g„ G M'^ such that 

\Ky-v^^'\\n + \qn\<2-\ (2.35) 
Let On = Y.k>n In: ^^^^ c„ and 

\\v^- - c„) - v-^\- - c„+i)||« < 2-^ (2.36) 

which implies t^^v^ — )■ v strongly in H, whence also f " — )■ w strongly in □ 
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We apply the mobile distance only to the continuous spectrum part because, on 
the one hand, the discrete spectral part is finite dimensional and smooth, and on 
the other hand, the linearized energy is conserved only on the continuous spectrum. 
Choose positive constants 6, Co, Ci, C2 such that 

< < K Co < Ci < C2. (2.37) 

The required smallness of C2S, I/Cq, Cq/Ci and C1/C2 is implicit in the following 
arguments, but only in terms of d, f, Q and k. Henceforth, we shall regard those Cj 
as being fixed constants and ignore the dependence on them unless it is important, 
while we regard 5 as a small parameter (with the smallness depending on C2), 
keeping track of its impact on the estimates. 

The quasi-distance : "H x "H — )• [0, 00) is defined by 

ms{v'',v'y := \\Pd{v'-v')\\l + ra^,{Py,Py)\ (2.38) 

where 4>s{ci) := (p{a/6) with a fixed (f) G C°°(M) satisfying 

M = [' O<0'<1. (2.39) 

la [a > 2C2) 

We will localize the equation for v within distance 0{S) from 0, such that the 
evolution outside of it becomes purely linearized. We have chosen 0^ such that the 
"fare" is purely proportional to the translation distance within the nonlinear region, 
but there is an additional fee for "excessive weight" over 0{6). It is easy to see that 
ffv^ has the same properties as m<^ in Proposition 12.21 

3. Construction of manifolds for a localized equation 

In this section we construct a global center-unstable manifold for a equation of 
V with localized nonlinearity around G "H. The manifold obeys the original fiow 
only on the subset 1-L± in a small neighborhood. 

3.1. Localization of the equation. Let x ^ C*^(M) be a non-negative symmetric 
decreasing function satisfying x(t) = 1 for |t| < 1 and x(^) = for \t\ > 2. Let 

Xs{v):=x{\\vrn/6'). (3.1) 

We will construct a center-unstable manifold near for the equation of v with the 
nonlinearity localized within 0{S) distance from 

vt = JCv + Msiv), Msiv):=X5iv)[A{v)-ViQ + v) + N{v)]. (3.2) 

Hence each component in the spectral decomposition solves 

dtXk± = ±kXk± + P±kMs{v), 

dtfi = -u + P^Ms{v), dtu = PuMs{v), (3.3) 
da = JC-f + P^Msiv). 
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Lemma 3.1. The equation fl3.2p is globally wellposed in %, and for any solution v, 

SMv\\v{t)\\E<\\vm\E. 

|t|<i ^ ^ 

sup|||7Wll|-||7(0)||||<(||t;(0)|UA5f. 
\t\<i 

Proof. Let f be a local solution around t = 0, and let w{t, x) = v(t, x — c), where c 
is the solution of 

c = X5{v)A{v), c(0) = 0. (3.5) 
Let Tc be the translation operator 

Tc(p{x) = ip{x - c), (3.6) 
then the equation for {w, c) is given by 

w = JVw + F{w,c), c = B{w,c), (3.7) 
with the nonlinear terms F and B, defined by 



where are translates of A, N: 



A,{w) := A(r» = {H{Q) - {VQ,\w,))-'uj{w,Vw)/2, 
N,iw) := r,iV(r» = /(Q, + w,) - f{Q,) - f'iQc)w,. 
Choosing some appropriate Strichartz norm, for example 



(3.9) 



|u.||s.,:=^ + (3-10) 



with p = -j^, we have by the Strichartz estimate for the free Klein-Gordon equation 

lkl||str{0,T) < ||w(0)||^+ ||F||l1W(0,T), \\c\\l°-{0,T) < |c(0)| + T\\B\\l^(^o^t), (3.11) 

where the nonlinear terms are estimated by Holder 
F(0,c) = 0, 5(0,c)=0, 

l<5K, c^)l < [hcKlklk A 5)2 + (||<t.1|« A 5m\w^\\n A 5), 

||<F(^/;^C^)|Ul^(0,T) < T[||<CnU^||t/7^||gtr + IHw^llstr] 

+ + ||<U»1|str](|k1|str A5), 

on the time interval < t < T <^ 1. Hence if 5, T > are small enough, we obtain a 
local solution of {w, c) on (0, T) in "H x by the standard iteration. By Gronwall, 
it is extended to any finite time intervals. In particular, 

lki||str(-i,i) < lk(0)lk, |c|l-(_i,i) < |c(0)| + i\\wiO)\\n A 5)^ (3.13) 
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and moreover, if ||ti;(0)||^ ^ 6 then ^ 6 for |t| < 1 and so, c(t) = 

and F = {0, f'{Q)wi). Hence we obtain by the usual iteration and continuation 
argument, 

ll<^1lstr{-l,l) < ||<^"(0)|k, 

hc^Lo^i-iA) < (ll<^'(0)||« A 6){\\w'mn A 5), (3.14) 

\\<F{w\cn\\LiH.i-i,i) <\\<w^mn. 

Next we prove the second and third estimates in (13. 4p . Since they are now obvious 
for ||f (O)ll^ > CqS, we may assume that 11^(0)11-^ < Co6. For the Vd part, we have 
by the energy inequahty 

\\vd-e'^%mL^E(o,T) < \\D-'Msiv)hin4,^T) < II^IIL(o,t) < HO)\\h- (3-15) 
For the 7 part, we have 

dtiCjh) = Xs{v)[A{v){f{Q)^,\V^,) + A{v){CVvdh) + {N{v)\^,)], (3.16) 
and so 

[il7lll]o < ||(A(t;)/'(g)7i, A(t;)D2t;,,iV(t;))|Uj^.(o,T)||7l|L,-W) 

<\\vmH 

and we are done. □ 
Denote the nonhnear propagator for equation (13. 2 p on H by 

U{t):n^n, U{t)v{0) =v{t). (3.18) 

3.2. Smallness of nonlinearity. The following estimate on the nonlinear term by 
the mobile distance will be the basis of all the succeeding arguments. 



Lemma 3.2. For any two solutions v^{t) = U{t)v^{0) E (^(M;?/) of (13. 2p . we have 

(3.19) 



supm5i;''(t) < msv''{0), 
\t\<i 

sup ||PrfK(t) - e-^^*<t;^(0)]|U +\[msPy]l\ < 6msv'^{0), 
l*l<i 



where the implicit constants are determined by d, f , Q and k. 

Proof. Without loss of generality, we may assume that for some q(t) G M.'^ 

rasv^itr = ||P.<t'^(t)|||+ ||/(t) - W(t)||| + |g(t)r0.(ll7^(t)IU)' (3.20) 
for \t\ < 1. Decompose each solution by 

= J2 K±9k± + li' -VQ + u^ ■ JVQ + -i^. (321) 

±,k&K 

(I) Case ||f^(0)||2 < Ci5: The previous lemma implies that ||f^(t)||2 ^ (^2^, and 
so 05(||7-'(t)||_E) = 1, for \t\ < 1 and j = 0, 1. The discrete components solve 

dt<f/ = -<z/^ + P^<Msiv''), dt<i^^ = P,<Ms{v^), ^ ■ ^ 
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where the nonhnear term is bounded by 

\\<Ms{v'')\\H-^<Sxhsv\ (3.23) 

which is proved as follows: (I2.29P and the translation invariance of u}{v, Vv) imply 
\<X»\<\<m\/S<msv''/6, 

\<Aiv^)\ < 5[\\D-'<v''\\h + ||7° - r.^'Wu] < Smsv'>, (3.24) 
h^V^Wn-- < \\D-'<in\n<xhsv\ 
The nonlinear part is estimated by using Sobolev 

\\<N{v'')\\h-^ < \\N{v') - T,N{v')h. + f \\deTe,N{v')\\^-.dd 



(3.25) 

< 5\\vl - TqvlWu^ + Iglll^^ii^i < 5m^^;^ 

where g := min(2, 1 + 1/p) for d > 2 and q := 1 for d = 1. Thus we obtain (13.231) . 
Therefore, we have for \t\ < 1, 

||<^>(t) - e'^^*<t;^(0)||E < f 5msv^{s)ds . (3.26) 
The linearized solution enjoys the obvious bound 



e 



''''<v''M\\E<e''\\<v'>M\\E. (3.27) 

For the difference in the 7 component, we need the mobile distance. Let 

r^=r,., C' :=rV, := r^Q = Q{- - c^), (3.28) 
for j = 0, 1, where c'it) G M.'^ are the solutions of 

= Xsiv^)Aiv^), c°(0)=0, c^(0) = g(0), (3.29) 
where q has been chosen in f l3.20p . Then we have 

C^' = jac' + ■ {P^^PK' - Pi^C') + ^^(0, ivp), (3.30) 

where O , Pj and P^^ are linear operators, and is the nonlinear part, defined by 
:= T^C{T^y, Pi := T^P^T^y, := xs{v^)N,j{z^), (3.31) 

and := r-'f^ + . Hence the difference satisfies 

<C = J/:°<C + P'{0, <M^) + R, \\R\\n < S[msv'^ + \<^\ + \\<Ch], (3.32) 

using fl3.24p as well as Hf-'H^ ^ 6. The Strichartz estimate for the free Klein-Gordon 
equation yields (regarding J(£° — 'D)<(^ as a perturbation, which can be done by 
partitioning the time-interval) 

ll<1lstr(o,i) < IK(0)||2 + ||<Mn|^i^.(o,i) + ||P||l-«.(o,i), (3.33) 

where the nonlinear part is estimated by applying Holder to the second order Taylor 
expansion of /: 

f{Q' + 4) - fm - rm^ = f !\nQ' + aez{)e{z{f] dade, (3.34) 

Jo Jo 
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and 

4f{Q- + z\) - fm - rmzi] 

[/"(g" + ez'^)<{Q'' + ez\) + /"(g° + aOzDoz^i^z^] dade, 

where := (1 — a)Q^ + aQ^ and 2;" := (1 — a)z^ + az^. Thus we get 

\\<M^LlLliO,l) < (1 + lk1|str(0,l))""'|k1|str(0,l)[||<C^||L-(0,l) + IH^II Str(0,l)] • (3-36) 

On the other hand, we have 

Ih^^'llstrCO,!) ~ IHc''||l°°(0,1) + ||<^dl|L°=W(0,l) + IHCII Str(0,l) , , , 

m5^'<ll<^^lk+IKIk + K|. ^ ■ ^ 

Combining these estimates with fl3.33p - fl3.36p . f l3.26p and fl3.27p . we obtain 

sup Xhsv" < ||<C^||ci(0,l) + ll<1lstr(0,l) + ||<^^dl|str(0,l) ^ mv'^iO). /g ggN 

For the sharper estimate on the 7 part, we use 
(m57^)2<||(r°)*<lll + Kr, 

\\{Ty<ai = {^'<\p^<n + \\Pdi\^ + <cnrE, 

with equahty at t = 0. For the distance term, we have from fl3.37p and fl3.38p 

|<c^(t)| < |<c^(0)| + CSmsv^iO). (3.40) 
For the translated part, we have 

\{c'<c\p^<c) - {c'<c\<)\ < \\<p':cx < l<c1^5^ 

\\P,^\x + <cn\\l<\<c-\'5\ 

and 

dt{C''<C\<C) = -{f'{Qco)c'' ■ VQ,o<K) + 2(£°<1^°<A^^ + R)- (3.42) 
Hence for |t| < 1, using fl3.38p and (13.360 as well, 

\[{c'<c\p^<cm 

Plugging this and f l3.40p into f l3.39p . we obtain the desired upper estimate on the 7 
part. For the lower estimate one reverses time, completing the proof in the case (I). 

(II) Case miuj ||f''(0)||2 > C06: The previous lemma implies that ||f''(t)||-H ^ 
C06 ^ 6, and so, for \t\ < 1 and j = 0, 1, 

v^it) = e'^'v\0), ||7^'(t)IU = ll7^'(0)|U, \\P,<v-mE<\\P,<v-mE. (3.44) 

To estimate the 7 part, let ((t) := 7°(t) — 'y^(t,x — q{0)). Then 

C - J£C = (0, [/'(Q) - f (Q,(o))]r,(o)7i), (3.45) 
where the right-hand side is bounded in by 

\qm\\l\0)\\n<6msv%0), (3.46) 
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where we used that a < 6(f)s{a). Hence using the energy inequahty for £, we get 

\m\OYo\<5xnsv^mC\\Lrn., 



llPdCWIU = WPdrml'llE < |g(0)|||7'(0)|k < Smsv^iO), 

and so 



(m.7^W)'<llCWIII+|g(0)|V.(liy(0)|U 



2 

< (m57"(0))' + C52(m5^>(0))2. 



(3.47) 



(3.48) 



This completes the proof in the case (II). 

(III) Case ||t;°(0)||^ > Ci6 > CqS > ||t;^(0)||^: The previous lemma implies 
that \\v%t)\\n > Ci6 > CqS > \\v\t)\\n for \t\ < 1, and so by ^TM . 

msv'-it) ~ \\v%t)\\H ^ \\v%0)\\n ^ xhsv'^iO) > 6. (3.49) 

For the difference from the linearized solution, (13.151) yields the desired estimate. 
The estimate on the increment of the 7 part is similar to the case (I), but now 
^0 _ ^0 gvolves linearly, which means that the nonlinear terms in <('^ depends only 
on (^^. Hence f l3.36p is replaced with 

IHM^IUji. = llMi^i^. < Wv'Wl, < 5\ (3.50) 

Noting that 5 < m5(f°(0), f ""^(O)), the rest of the argument goes through as in case (I) 
above. 

(IV) Case ||w°(0)|| < Cq5 < Ci5 < \\v^{<d)\\u. Although this is symmetric 
with the previous case, we have to check the mobile distance part, since there we 
introduced asymmetry with (13.201) . The difference appears in (I3.39P : 

{mv-f < ||(r°)*<1ll + \<e?MW\\E?. (3.51) 

However this is admissible, because 

[\<c'\UH\\e% < 6'M\\v\0)\\e) ^ Smsv>{0), (3.52) 

and the remaining argument is the same as in the previous case. □ 

3.3. Evolution of graphs, center-stable case. Now we consider the graphs of 
f>o ^ f _ satisfying a Lipschitz condition. It is convenient to extend them to the 
whole H. Our class of graphs for the contraction argument is given by 

ge,s ■= {G -.H^ P-H I G = G o P>o, G'(O) = 0, ||<G(t;^)|U < ^m^^^^} (3.53) 
for small £ > 0, and the graph of G G Q^^s is denoted by 

\G\ := G -H I P^v = G(v9)}. (3.54) 

A center-unstable manifold will be found as the unique invariant graph, by the 
contraction mapping principle in Qi^s. 

For p > d/{d — 2), the Sobolev inequality does not imply that v is bounded in 
L^, and consequently we can not prove strict invariance of Qg^s for t > 0, but the 
"almost invariance" given below is sufficient for the contraction argument. 
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Lemma 3.3. There exists Ci>l such that ifi,6>0 satisfy 

W + 6<t:k, 5 < ek, (3.55) 
then for any two solutions v^{t) = U{t)v^{0) (j = 0,1) satisfying 

\\<v''-iO)\\E<ixhsv''{0), (3.56) 

one has 

fe:;:;"' irj^;:,,. 

Proof. The linearized solutions of the discrete modes are estimated by 



min(e±^*,e±^*)||P±(/.|U < Il^^ie-^^VIU < max(e±^*, e±'=*)||P±(/.|U, 



e 



-2«t(l _ £2) ^ g-2H^2 _ C6^]^SV^(^Q)2 (Q < t < 1), 



(3.58) 



(3.59) 



The previous lemma implies that 

\\<v^-{t)\\E < ||<e^^V(0)||E + C5m5^^(0) 

< [max(e-^*,e-^*)£ + C5]m5W^(0), 
and also, 

mv^ity > ||e-'^*<T;^(0)||| + m,7^(0)2 - C(5W(0)^ (3.60) 
Plugging fl3.58p into the last estimate, we obtain 



*^""w^£{;.-..„l_p,^,^^,„.(„,. (3.61) 

Combining it with f l3.59p yields for \t\ < 1, 

||<^-(t)||E < (e^£ + C5)(e^ + C5)m5t;^(t) < 2ie^^msv^{t), (3.62) 
provided that 5 < £ < 1. For 1/2 < t < 1, we obtain 

\\<v''-{t)\\E < (e-^/2^ + CT)(e-2«^i _ ^ ^-2kf _ c6Y^/^msv^ 
< (1 - k/3 + C6/i)i[l + C{6 + K + fk)]msv'' < imsv", 
under the condition fl3.55p and k, <^ k < 1. □ 



(3.63) 



As an immediate consequence of the above lemma together with a mapping degree 
argument, we obtain the following result. 

Lemma 3.4. Under the condition f l3.55p . U{t) for \t\ < 1 defines a map U(t) : 
Qi,5 — ^ GclI^^ uniquely by the relation U(t) \G\ = \V{(t)G\ . Moreover, if 1/2 <t < 1, 
then U{t) maps Qi^s into itself. 
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Proof. The previous lemma yields for any ^9°, ip^ E U {t) \G\ , 

IHV'-IU < (3.64) 

Since £ < 1, it implies I^v^^Hb < m^v?^- Then the conditions U{t)\G\ C \Kit)G\ 
and U(t)G o P>o = U{t)G define U{t)G uniquely and consistently on the set 

P>of/(t) \G\ + P^V.. (3.65) 

The proof is complete once the above is shown to be for which we use the degree 
argument. Suppose for contradiction that there exists ijj G P>qH \ P>oU(t) \G\. In 
other words, for any a G P-Ti, 

U{-t){a + ij) ^\G\. (3.66) 

Let m(a) := P_U {—t){a + — G{U {—t){a + tp)), then m is a continuous map from 
P-H to itself, such that ^ m{P_'H). On the other hand, if \a\ ^ 5, then 

m(a) = e-^'^'a - G'(e-'^^V)- (3.67) 

Define $ : ^ P_H, ^ : \ {0} ^ 5^'-^ and m : (0, oo) x S^-^ S^-^ by 

$(X) = ^Xfc^fc_, vI/(X) = ^, m(P,^) = v[/o$-iomo$(P^). (ggg) 

Then m is continuous, but the degree of m(i?, ■) is for small i? > and 1 for 
large i?, which is a contradiction. Hence U{t)G is well-defined as a map on "H which 
is right-invariant for P>o. The Lipschitz bounds are immediate from the previous 
lemma. □ 

3.4. Contraction of graphs, center-stable case. We introduce the following 
norm in ^ := {jf^^^Qi^s 

G g .= sup . (3.69) 

It is easy to see that the set Q is independent of 5 > 0. For any G G Q^^s and any 
G "H, we have ||G'(V^)||ij < £m5(V^,0) < ^HV^IU, and so 

\\G\\g < L (3.70) 

^ is a Banach space with this norm, where each Qn^i, is a bounded closed set. The 
contraction argument is completed by 

Lemma 3.5. In addition to fl3.55p . let 

6<^k\ (3.71) 

Then the map U{t) is a contraction on s for all t > 1/2. 

Proof Let T G [1/2, 1]. For any G^ G Gi^s for J = 0, 1 and any eU, let 

y^ (t) := Uit - T)[P>o^ + (W(T)G^>]. (3.72) 

Since P>of°(T) = P>ov^(T), we have 

m5t;^(T) = ||<t;^_(T)|U. (3.73) 
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Applying Lemma [3.21 from t = T, we get for < t < T, 

||<^>(t) _e^A*-T)^^> + < S\\<v'>_{T)\\e. (3.74) 

Hence using the same estimate on the hnearized solution as in fl3.58p 

hv^^miE > ||e~-^^W (r)|U - C6\\<v^_{T)\\e > {e^ - C6)\\<v^_{T)U . . 

On the other hand, since vi(0) = G^{y^{Q)) and G ^^,5, 

ii<t^'-(o)iu < hG^{v\m\E + hG\v-{m\E 



< hG^\\g\\vlM\\E + ^^8vlM- ^^'^^^ 
dM]) as well as (1338|) yields 

ll4o(0)ll|<(e^^^ + C<5)||^|||. (3.77) 

Inserting this and the second inequality of fl3.75p into fl3.76p . we obtain 

hv^-mE < {e^^ + GV5)\\<G^\\g\mE + m<v''.{T)\\E. (3.78) 
Combining this and (13. 75 p . we conclude that 

\\<v'-_{T)\\e < (1 - G5i)-\e^^ - G5)-\e^^ + GV~5)\\<G^\\gU\\E 



(3.79) 

<e-(^-«)^(l + C%/5)||<G^||a||V'|U. 



f l3.7ip and n <^ k imply that there exists A < 1, determined by A;, 5, £ such that 

< AING^II.. (3.80) 

\mE 

Taking the supremum over alXip El-L yields 

||<W(r)Gnig < AIHG^Ie, (3.81) 
as desired. The case T > 1 is now obvious by iteration. □ 
Thus we obtain 

Theorem 3.6. Suppose that £, 5 > satisfy fl3.55p and f l3.7ip . Then there exists a 
unique G Qi^s such that U{t)G^ = for allt>0. The uniqueness holds for any 
fixed t > 0. 

Proof. For any T > 1/2, the above lemma implies that there is a unique fixed point 
of Z//(T) in Qi^s- Since the equation is invariant for time translation, it implies that 
U{t)G G GcLe,s is also a fixed point for all < t < 1. Then the uniqueness of the 
fixed point implies that U{t)G = G for all < t < 1, and so for all t > 0. If 
U{t)H = H for some t > and some H G Qe,^s, then by iteration U{T)H = H for 
some T > 1/2, and so H = G. □ 
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Since U{t) is invertible, U{t) \G^\ = \G^\ for all t e R. 
The conditions fl3.55p and (13 .71 1) are satisfied for £ = 0(5) as 5 — )■ +0, which implies 
that 

\G,\ 3 ip, Mn < S =^ H^,gk+)\ < 5\ (3.82) 

in other words, [G*] is normal at to (—A;, ^)pk for each k E K. 

Notice that the above construction did not really use the special property of the 
generalized null space of the linearized operator. However, the constructed manifold 
makes sense for the original equation only on the subset 1-L±, for which we need the 
property that the generalized null space is exactly generated by the symmetries of 
the equation. 

3.5. Evolution of graphs, unstable case. We now carry out an analogous proce- 
dure for the finite-dimensional unstable manifold. Thus, we now consider the graphs 
of f+ I—)- f <o satisfying a Lipschitz condition 

g+s := {G : -H -> P<o'H I G = G o P+, G(0) = 0, msG{v^) < £||<<|U} (3.83) 

for small £ > 0, and the graph of G G Q'^g is denoted by 

[GJ ■.= {ven\ P<oV = G((/.)}. (3.84) 

The unstable manifold will be found as the unique invariant graph, by the contrac- 
tion mapping principle in Q^^. We formulate the analogue of Lemma 13.31 in this 
case. 

Lemma 3.7. There exists Cl > 1 such that if £,S > satisfy (I3.55p . then for any 
two solutions v^{t) = U(t)v^{0) (j = 0,1) satisfying 

m5<o(0)<^||<<(0)|U, (3.85) 

one has 

A v" < /Ci^lN^WIlE (l«l<i). rtsfi^ 

- \t\\..-4mE (i/2<i<i). 

Proof. We again have (I3.58p for the linearized solutions of the discrete modes. In 
particular, 

||<<(t)IU > min(e^*,e^*)|K(0)b - C5rasv'i0) (3.87) 



(3.^ 



By Lemma [3. 2[ for t > 0, 

msv^<o{ty < he'^-X^ml + m57"(0)' + G5W(0)' 

< {e^'^'e^ + G6\l + e^) + ^')||<<(0)||| 

and one now concludes by combining these estimates, cf. Lemma [3.31 □ 

One now has the following analogue of Lemma 13.41 

Lemma 3.8. Under the condition (I3.55P . U{t) for \t\ < 1 defines a map U{t) : 
Qts ~^ ^CLiS uniquely by the relation U{t)\G\ = \U{t)G\. Moreover, if t > |, then 
lA{t) maps into itself. 
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Proof. The mapping properties for |t| < 1 and ^ < t < 1 are an immediate con- 
sequence of the previous lemma. The extension to t > | then follows by itera- 
tion. As in the case of the center-stable version, the main issue is to show that 
P+U{t) \G\ = P+n = V.+ for all \t\ < 1. Thus take e V.+ and denote the R-hsl\ 
in Hj^ by B^. Lemma [3.21 implies that if < £ ^ 1 then 

||p+f/(t)(^ + emu >R V |t| < 1, G dB+ 

for any G G Q^^, and with absolute implicit constants. This shows that, with 

deg(P+[/(t)$,5+,V'o) = l V|t|<l 
provided R is sufficiently large, and we are done. □ 

3.6. Contraction of graphs, unstable case. Let := [Ji^qQ^s- before, the 
set is independent of 5 > 0. We introduce the following quasi-distance d+ in Q~^: 
for any G\ G 0+ let 

d+iG\G') ■.= snp '^fy\ (3.89) 



It is clear that this expression is finite, and that it satisfies a triangle inequality with 
the same multiplicative loss as in Proposition 12.21 

d+{G\ G^) < Cd{d+{G\ G^) + d+(G^ G^)). (3.90) 

Moreover, is a complete quasi-distance space, in which Q^^ is closed. Recall that 
the Banach fixed point theorem is valid in complete quasi-distance spaces: 

Lemma 3.9. Let X be a complete quasi-distance space, and let A : X ^ X be a 
contraction. Then there is a unique fixed point E X of A, which is obtained by 
x^: = lim„_>.oo ^"(a;) for any x E X. 

Proof. Let G > 1 be the constant in the quasi-triangle inequality in X, let A G (0, 1) 
be the Lipschitz constant of A, and fix m G N so that A™G < 1. Take any xq G X 
and let Xn = A"™(xo) for each n G N. Then 

d{xn+i,Xn) = rf(A"^(x„), A"^(x„_i)) < A™d(a;„, Xn-i) <■■■< A"^"d(Xi, Xo) (3.91) 
Hence for any k > j > 1, hj repeated use of the quasi-triangle inequality, 



d{xk,Xj) < Cd{xk,Xj+i) + Cd{xj+i,Xj) < ■ ■ ■ < G' ^d{xi,xi^i] 

i=j+i 

h 

< (GA™)'"^A-^rf(xi,Xo) < ^— ^d(a;i,xo). 
i=j+i 



(3.92) 



Hence Xn 3xoo G X, and by the continuity of A, A™(xoo) = Xoo- Then 

rf(A(Xoo),Xoo) = rf(A'"+l(Xoo), A'"(Xoo)) < A"rf(A(Xoo),Xoo), (3.93) 

which implies that A{xoo) = a^oo- The uniqueness follows in the well-known way. □ 
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The following is an analogue of Lemma I3.5[ but here the evolution time has to be 
long enough to absorb the quasi-triangle factor in using the "chain-rule" in Q^. 

Lemma 3.10. There are So > andT > such that if 6 < 6o and f l3.55p is satisfied, 
then the map U{t) is a contraction on Q'^g for all t > T. 

Proof Let G\ G 0+^, T>0,iljeH+, and 

v^(t) := Uit -T)[i; + (W(T)G^>], (j = 0, 1). (3.94) 

Iterating Lemmas 13 . 1 1 and 13.21 from t = T down to t = 0, we obtain 

< Ce^^llV^IU, mv'-it) < Ce^^m5t;^(r) = Ce^'^xasv^iT), (3.95) 

for < t < T with some constant C > 1 (determined by d, Q, f and k). Hence if 

6 « e"2^^/C2, (3.96) 

then by iteration of those lemmas again, we deduce that 

K(o)IU<e~^^/'II^IU, 



(3.97) 



£5' 



(3.98) 



Since <o(0) = G^{v^{0)) and G^ E 

m5<o(0) < G4vasG'-iv\0)) + vasG\v''iO))] 

<G4d^iGn\\vlmE + i\\<vliO)h 
Plugging f l3.97p into the above, we obtain 

rhsv%o{T) < e'^^G,[d4G'')Ge~!^^/^^\\E + i^/Smsv%{T)]. (3.99) 
Choosing T so large while keeping f l3.96p . we can ensure that 

m,(W(T)G^)V^ = msv%{T) < AH^-IU, (3.100) 

for some constant A G (0,1) determined by d, f,Q, k,T and 6. Obviously, this 
remains to be true even if we replace T with any T' G [T, 2T] , taking 6 even smaller 
if necessary. Hence iterating V({t) allows one to draw the same conclusion for all 
t>T. D 

By the same arguments as in Theorem 13.61 one now concludes the following. 

Theorem 3.11. Suppose that i,6 > satisfy the assumptions of the previous lemma. 
Then there exists a unique G^ G Q^g such that U(t)G^ = G^ for all t > 0. The 
uniqueness holds for any fixed t > 0. Moreover, ifv(fS) G \G^\, then ||f/(t)w(0)||-^ — j- 
exponentially as t ^ -co; in fact, for any e > 0, 

e-'^^-'^'\\U{t)v{0)\\n ^ t -oo (3.101) 
Proof. The estimate ( I3.10ip follows from the previous proof. In fact, (13.970 implies 

IKi)IU<e5*||t;(0)|U (t-^oo), (3.102) 

where v{t) := U{t)v{0), but we could take the exponent arbitrarily close to k by 
choosing S even smaller. Since v{t) comes into any 6 ball as t — —oo, we may apply 
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such decay estimates for t sufficiently close to — oo, thereby deducing (13.1011) for all 
e>0. □ 

4. The unstable manifold 

We now describe the unstable manifold in the original u-formulation of the equa- 
tion, see fll.ip . Let be as in Theorem I3.11[ For any f (0) G with 

\\v{0)\\n < 6, define v{t) = U{t)v{0), 

uit) = {Q + v)it,--cit)), c(0) = co, cit) = A{vit)) (4.1) 

where cq G M"' is a fixed vector. By construction, u solves (II. ip . and by (IS.lOip one 
has c(t) —7- and c{t) — )■ c(— oo) exponentially fast as t — j- —oo. In particular, u has 
vanishing momentum: P{u) = 0. Then by design (cf. (I2.24p -( l2.26p ). uj{v, JVQ) is 
constant, and since it converges to zero as t — )■ — oo, must vanish. To summarize, 
we have obtained the following characterization of the unstable manifold. 

Corollary 4.1. The unstable manifold M.u is the set of all 'u(O) with u defined in 
terms of by means of (14. ip . Ai^ is invariant in backward time, and all solutions 
starting in M.^ converge to a trajectory of the form Q{- — c{t)) exponentially fast as 
t — 7- — oo, with c(t) —7-0 as t ^ — oo exponentially fast. J^u is a Lipschitz manifold 
of dimension K + d. 

The dimension count is a result of the fact that [G^J is of dimension K, and the 
translations (see Cq in (14.11) ) add another d dimensions. 

5. Trapping property of the center-stable manifold 

5.1. Restriction by the orthogonality. For any Banach space X, denote the ball 
around of radius i? > by 

BR{X):={^eX\M<R}. (5.1) 

Lemma 5.1. If i < 1, and 6 > is small enough (depending only on d,f,Q), then 
for any G G Qi^s, there is a unique map G : P.y+BsiTi) — > Pq^'H such that 

\G\: = {^p + Gi^)\^eP,Mn)} 

= {ipE \G\ n -Hx I P^+v e P^+Bsin), \PM <6}. 

Moreover, G is Lipschitz continuous in the mobile distance xhs. 

Proof. For any ip and z/ G M"^, put 

V^(i/) := + z/ ■ JVg, ^{v):=ij{u) + G{ij{u)). (5.3) 

It suffices to show that for any ip G P^^BsiTi), there is a unique fixed point u G 
Bs(R^) for the map 

u ^ M{v) := H{Q)-^uj{ip{u), V^(z/))/2. (5.4) 
For any G BsiW^) (j = 0, 1), we have 

<^(zy>) = <zy> . jvg + <G{iIj + ■ JVQ). (5.5) 
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Hence using WipWu < 5 and G G Qi^s as well, we deduce that 

W{^)\ < WmWl < 5^ l<^(^(^^), V^(z/^))| < 5\<u^\. (5.6) 

Therefore M : B^iW^) — )■ B^i^'^) is a contraction for small 5, ^ > 0, and so has a 
unique fixed point v = M^u) G i?5(M'^). Since M is Lipschitz for ip in the mobile 
distance, so is the fixed point ^{ip)., as well as G. □ 

Hence \G\ fl "Hj, is a Lipschitz manifold in the mobile distance around with 
CO dimension K + 2d = dim Po^T-L. 

5.2. Solutions on the center-stable manifold with the orthogonality. Let 

G = G* G Qe^s be the map for the center-unstable manifold of the localized equation 
given by Theorem 13. 6[ and let G = be the map for its orthogonal restriction 
given by the above lemma. The invariance of \G\ means that for any v{0) G [GJ, 
v{t) := U(t)v{0) stays on [GJ for all t G M. Let c(t) be the solution of 

c(0) = 0, c(t) = A{v{t)), (5.7) 

and u(t) = {Q + v)(t, x — c(t)). If f (0) G Bs{'H), then u solves the original equation 
(11. ip as long as v{t) G BsiTi). Meanwhile, the momentum P{u) = a;(f,VQ) + 
ijj{v, Vf)/2 is preserved, and so is co{v, JVQ), because of c = A{v). Hence if f (0) G 
Ti^ n BsCH), then v(t) remains there as long as v(t) G BsiTi). 

To see that the solution stays in the neighborhood for t < 0, expand the conserved 
energy by m = {Q + v){x — c) 

E{u) - J{Q) = -J2 kh+Xk- + liCill) - C{v), (5.8) 

where the nonlinear energy G is defined by 

Giv) := f{Q + v^) - f{Q) - ifrnv,) - ^{f"{Q)v,\v,) = o(||t;||^). (5.9) 
Suppose that for some to < 

\\v{0)\\n = e <^ 6, maxJv{t)\\H < 6, \\v{to)\\n :^ Vi6 + e{<^ 6). (5.IO) 

Then \E{u) - J{Q) \ ~ e^. Since v{t) G [GJ and G G Qe^s, (EID implies that 

mi) < E ^1 VA.-I + E{u) - J{Q) < £6' + « 5^ (5 ^) 

k 

for to < t < 0, which together with v{t) G 'H± implies 

^|A,+ (to)P^||t;(to)|lH- (5.12) 



k 



Now consider the nonlinear energy functional 

2' 



£{v) : = E{u) - J{Q) + J2 |( V + Aa 



keK 



(5.13) 



2 
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Since v(t) G Ti^, we have £{v) ~ \\v\\y^, and moreover, the equation of Xk±, see (I3.3p 
together with conservation of E{u) yields 

j^S{v) = J2 - At) + o{\\v\\l\\,^\), (5.14) 

k 

and so 

j^S{vm>Y,k'\l^^S{v{t,)). (5.15) 

k 

Therefore £{v{t)) cannot increase beyond 0{15^ + e"^) as t < decreases. 

In conclusion, for any f(0) G |"G*J such that ||f(0)||-H ^ 5, the solution v{t) = 
U{t)v{0) remains on [G*J, ||t;(t)||w < ^ for alH < and u{t) = {Q + v){x - c) 
with c = A[v) solves the original equation (11. ip for all t < 0. Thus we obtain a 
center-unstable manifold of the original equation with zero total momentum. 

More precisely, fix < 5' ^ 5 and let 

A^cn,o := {(Q + V{t)^){x - c) I G [aj , W^pWn < y, t < 0, c G R'^}, (5.16) 

then for any initial data 'u(O) = (Q + U{t)(p){x — c) G A^cu,0; the solution u(t) of 
flTT]) is on Mcu,o for alU < and P(M(t)) = 0. Moreover, U{t)(p G [G*J and 
||t/(t)(/9||-^ < 5' for all t <0. The nonlinear projection 

P±:no3 u^v en±; u = {Q + v){x-c) (5.17) 

is uniquely defined in a neighborhood of the translation family of stationary solutions 

MQ) ■■= {Qi^ - q)},m^ C -Ho, (5.18) 

by solving the equation 

= uj{v, JVQ) = uj{u{x + c), JVQ) = uj{u, JVQix - c)). (5.19) 
Indeed it can be solved locally by the implicit function theorem, since if u = Q{x — 

VMu, JVQ{x - c)) = -uj{Q{x - Co) + V, JV^Qix - c)) 

= H{Q) + Oi\c-co\+6). 

Since the mapping u ^ c thereby defined is smooth, the map mi—)- (f , c) is (locally) 
bi-Lipschitz in the mobile distance from T-Lq to 'Hj_ © M.'^. Since A^cm,o is mapped 
onto a 0-neighborhood of [G*J © M*^, the codimension of Aicu,o in is equal to 
that of [G*J in "H^, which is K. 

Thus we have obtained a center-unstable manifold A^cm.o of ^o(Q) in with 
codimension K. Its time inversion 



-Mcn.O =: Mes,0 C (5.21) 

is a center-stable manifold of J^oiQ)- 
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5.3. Lorentz extension of the center-stable manifold. Using the Lorentz trans- 
form 

u{t,x) y-^ Up := u{{p)t + p ■ x,x + p{{p) - l)\p\~'^p ■ X + tp) (p G M"*), (5.22) 

we can further extend A^cs,o to a manifold M.cs of codimension around the soliton 
manifold y{Q). Indeed, (11. ip is invariant for any Lorentz transform, while the total 
energy and momentum are transformed 

E{up) = E{u) (p) + P{u) ■ p, P{up) = P{u) (p) + E{u)p. (5.23) 

E"^ — |Pp is invariant, which is positive around S^{Q). Hence there is a unique 
p G M'^ for each solution u near the traveling waves, such that P{up) = and 
E{up) = ^Eiuy-\PiuW. 

However, one needs to be more careful because the Lorentz transform mixes space- 
time and a solution from Mcs,o may not be global in the negative time. Indeed, 
from [TTl [T2] we know that "half" of the solutions on A4cs,o (namely, as given by the 
separating surface Aicufi) blow up in negative time, at least when Q is the ground 
state and f{u) = p> 1 + A/d. 

The local wellposedness implies that for any T > there is 5 > such that for 
any initial data within distance 6 from ^o{Q), the solution extends at least for times 
|t| < T. The exponential decay of Q implies that for any S > there is i? > such 
that for any such initial data, the free energy in the exterior region |x — g| > -R is 
less than 0{6^) for some q G M'^. The local wellposedness, the conservation of the 
energy and the Sobolev inequality implies that every solution with small initial free 
energy is global, keeping the same size of free energy for all time. Hence the finite 
speed of propagation of the free Klein-Gordon equation implies that for small 6 > 0, 
every solution with free energy 0{6^) in |x — g| > i? at t = is extended to the 
whole exterior cone |a; — g| > R+ \t\ with the same size of the free energy on any 
time slice of it. 

Thus in conclusion, there is i? > and S{T) > for any T > such that every 
solution starting on Aics,o and within distance 6 from ^o{Q) is extended to the 
space-time region 

{{t,x) G M^+'^ I t > -T or \x-q\ >i?+|t|}, (5.24) 

for some g G M'^. For any Lorentz transform L, there is T > such that the image 
of the above region under L contains {{t,x) \ t > 0}. In other words, the image of 
any solution on A^cs,o close to ^o{Q) is extended to a forward global solution. The 
invariance of the solution set of Mcs,o for the space and backward time translations 
is also inherited by the image, because such a translation of the Lorentz transform 
is the Lorentz transform of another translation. It is also easy to see that these 
solution remains close to the corresponding traveling wave. 

However, it seems difficult to make the above argument uniform with respect to 
the Lorentz transform: the larger the momentum p, the smaller the neighborhood of 
^o{Q) needs to be chosen. This is why the resulting manifold is not strictly Lorentz 
invariant, but only within a neighborhood of ^{Q) depending on the Lorentz trans- 
form (but the neighborhood can be chosen uniformly for p in compact sets). 



28 



K. NAKANISHI AND W. SCHLAG 



Thus we obtain a center-stable manifold Aics of the soliton manifold ^{Q). Aics 
can be identified with the set of forward global solutions starting from it, where each 
solution is characterized uniquely by the total momentum and its Lorentz transform 
with 0-momentum starting from Aics,o- In this wajH, we can define a bi-Lipschitz 
map from a neighborhood of ^{Q) in H to a neighborhood of ^o{Q) © in 
'Ho©IR"'- Since it maps Aics onto the intersection of Aics,o®^'^ with a neighborhood 
of yoiQ) © the codimension of Mcs in "H is also K. 

5.4. Solutions off the center-stable manifold. It remains to describe the dy- 
namics off the manifold, or more specifically, the repulsive property of the center- 
unstable manifold in negative time. For this we need some sort of opposite to Lemma 

Lemma 5.2. Ifi,5>0 satisfy 

S < ik, (5.25) 
then for any two solutions = U{t)v^{0) (j = 0,1) satisfying 

max(||<t;g+(0)|U,m57"(0)) < £||<t;"_(0)|U, (5.26) 

one has 



and 



max(||<.>,(t)|U,m,7^(t)) < { J"'/' < ' " °K (5.27) 

VII o+v;iii?, M V _ ^^ii^^.^^^ii^ (-1 < t < -1/2), ^ ' 



\\<v^ mil. > J (-1/2 < t < 0), 

"^'-^'^"^ - \e-^*/^||<^;^_(0)IU (-1 < t < -1/2). ^^^ ^^ 

Proof. Let m(0) := ||<t''5l(0) H^; ~ msv^(0). Lemma |3T2] implies that 

\\<vl4t)\\E < (e-"*£ + C5)m(0), \\<v''_{t)\\E > (e"^* - CT)m(O), 
ms^^it) < (£ + C(5)m(0), 
for — 1 < t < 0. Hence 

max(||<t;^+(t)||s,m57"(t)) < (e'^^f + C<5)(e-^* - C5)-i<t;^_(t)|U 

<e(^-)*(£ + C(5)||<t;^_(t)|U, 
and the conclusion follows from (15.251) as well as k <^ fc. □ 
Let w°(0) G Bs{n) nn± \ \G,\, and 

i> := P^+v%0), v\0) ■.= i> + G,{ij). (5.31) 

Then we have 

<<+(0) = 0, ||X(0)||| < ||<.;^(0)||| ^ hviml < 8\ (5.32) 



(5.29) 



(5.30) 



^The general solution close to S^{Q) may well blow up in both time directions, but the smaller 
neighborhood yields the bigger lower bound on the existence time, which is sufficient for the 
construction of this bi-Lipschitz map. 
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Hence we can repeatedly apply the above lemma to deduce that 

\\<vimE>le-!^'/^<v'LmE (5.33) 

for all t < 0. In particular, 

wv^mE > ii<^-(t)iu - h-imE » 6 (5.34) 

for sufficiently large —t. 

In short, any solution starting from "H^ fl BsiTi) \ [G*] moves out of the neigh- 
borhood Bs{l-i) for large — t. Of course, this is meaningful for the original equation 
only until the (backward) exiting time, but it implies that any trapped solution in 
'H±_ within distance 5 must be on the manifold [G^,J for large —t. 

Combining this with the result in the previous section, we conclude that the local 
center- unstable manifold A^cu.o is characterized as the collection of solutions with 
0-momentum which stay close to ^o(Q) for all ~^ ^ 0. By symmetry, M.cs,o is the 
collection of solutions with 0-momentum which stay close to S^oiQ) for all t > 0. 

Let A^cu,o be the maximal forward evolution of A^cm.O; and let A^cs,o be the 
maximal backward evolution of M.cs,o- Then M.csfl is the collection of solutions that 
stay close to ^oiQ) for large t, namely the initial data set for which the solution will 
be trapped by ^o(Q)- We have the same characterization for A^cu,o for t — oo. 
By the Lorentz transform, we can extend them to solutions with nonzero momentum 
which are trapped by S^{Q) with the same momentum. 

6. Regularity of the center-stable manifold 

The above construction implies only Lipschitz continuity of the manifold. For the 
differential structure of \G^\ , we also have to take account of the spatial translation. 
In the following, we assume that / satisfies f ll.7p and a := max(l,p — 2). 

Definition 6.1. Let y be a Banach space. We say that a function G : H — t- F is 
mobile- dijferentiable at G 7/, if there is a bounded linear M : x R'^ ^ Y such 
that 

lim||G(^(,)) - [Giip')+sMi^,q)]\\/e = 0, (6.1) 

where ip(^^) := {ip + eip){x + eg), for any (■?/', g) E Ti x M'^. It is obvious that M is 
unique. We call &G{(p) := M the mobile derivative of G at ip. 

Let G'{(f) be the usual derivative in the Frechet sense. Then we have 

^G(v9)(^, q) = G'(v9)(^ + ■ q), (6.2) 

provided that G is differentiable in the VH topology, but in general, it makes sense 
only in the subspace q = 0. Hence the mobile-differentiability is stronger than the 
differentiability in "H, and weaker than that in VH. 
If G G Qe,s and mobile-differentiable, then 

||G(^(,)) - G{if')\\n < ira^im^ v') < HMn + \q\MMn)], (6.3) 
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which imphes 

\\^G{^){^, q)\W< + WsiMW)]- (6.4) 

Moreover, we have 

G{{ip + e^ + o{e)){x + eq + o{e)) = G(<^(,)) + o{e) 

= G{^)+e^G{^){ij,q) + o{e). ^ ' ^ 

We are going to prove that G^. : l-L P-T-L is "mobile-C^'°", by showing the 
flow-invariance of the following set of such graphs. 

Definition 6.2. For each 6,i,A > 0, and a G (0, 1], we define Q"^'^ as the set of all 
G G Qi^s that are mobile differentiable at every (p ETi, satisfying 

||^G(y.°)(V^,g)-i^G(^i)(r*^,g)|U 



for all (f^jip ^'H and q,b E M'^, where 



x"'' := Ixl"" + \x\. (6.7) 



We will prove that Q"^'^ is invariant by the flow, provided that 6, £ are small and 
A is large. First we investigate the backward evolution of the mobile derivative. 
Assuming the smallness of ^, 5 > as in ( I3.55P and f l3.7ip , for any G G Qi^s and 
t > 0, define Gt : H ^ P-U and Gt : U ^ U by 

Gt:=U{t)G, Gt{v)--=V>^ + Gt{'^). (6.8) 

Let iljen^q^be R'^, and to e [0, 1]. For small e G M, let 

V{s){t) ■=U{t-to)Gt^-,{ip^e)), W(^e){t,x) =V(^e){t,X-C^e)), (6.9) 

where (w,c) = (w(£),C(e)) is the solution of (13. 7p with the initial data 

W(e){to) = {f + eip){x-b), c^^^{to) = b + eq. (6.10) 

Since the nonlinear term {F, B){w, c) in (13. 7p is G^ from Str x L'^ to LjH x Lf^, with 
a small factor on a short time interval (0,T), it is straightforward by the iteration 
argument that {w, c) is differentiable in Str x at e = 0, with the derivative 

lk||str(0,l) + IklU-CO.l) < lk(to)||K + lfi'(^o)|- 



(6.12) 



Let 1] := z{t, X + C(o)(t)) and F = F — {0, f'{Qc^^^)wi). Then 

7] = JCri + (0, g ■ Vf{Q)vr) + B{w, c) ■ Vr^ + r;9(^,,)F(«;, c) ■ (2;, g), 
r]{to) = P>o(^ + g ■ V<^) - g ■ VGtM + S^GtMii^^ 9{to) = q, 

where the subscript (0) is omitted. Mobile-differentiating the identities 

P„%)(to) = Gt^iifie)), P-V(e)it) = GtiVi,)ito)) (6.13) 
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yields 

PAvito) + q ■ VG„(^)] = !^GtMi^, l), 
P_[r/(t) + git) ■ Vvit)] = &Gt{v{t)){r^{t),g{t)). 

Since G{ip) = G{P>oip) and 

P>o^{e) = P>oKv>o + ^^>o)(a; + eq) + eq ■ VipJ\ + 0(e), 

we have 

^G(<^)(^, q) = ^G(<^>o)(^>o + q ■ P>oV^., q). 

Lemma 6.3. Let £,5, A > satisfy ( 1X53]) . (1X7T]) and 

A > 1/5. 

Then for any G G Q'^'^ and any t^ G [1/2, 1], G G^'^ ■ 

Proof. First of all, ( 16. 4p is enough to have (16. 6p in the case where 

£/A« ||<^o-r,(^i||^ + |6|0,(<^i). 

Hence we may assume 

m^<^^ < 1|<^° - wii^ + |6|</>5(/) < ^/A « 5. 

Therefore we have either ||v5°||^ — ^ ^ or ||v5°||-h + llv^"^l|w ^ ^■ 

Next we investigate evolution of the difference of mobile- derivatives. 
(^°, <^^,ip en and q,be R'^, let 

t;°,)(to) := G,„(r;(/ + e^)), t;(\)(to) := G,,(r,;(^i + er*^)), 

for j = 0, 1, where c|^^ is the solution of 

^(e) = Pi^U^ 4))' 4)(^o) = ^g, 4) (to) = b + eq. 
Let {z^,g^) be the derivative at e = of (w;^^^, c^^-j): 

^^• = lim4^. 
Henceforth the subscript (0) will be omitted. The initial values are 

w\to) = Gt,{y:>'), c°(to) = 0, w\to) = TAo{y^'), c\to) = b, 
z\to) = P>o^ + q ■ Gto(<^°) + ^Gto(<^°)(V^, <l), 9\to) = q, 
z\to) = + n[q ■ Gt,{v') + ^Gi,(^i)(r*^, g)], g\to) = q, 

where P^ := nP^r^, and Gt : "H "H*^ is defined by 

= VP_y^ - P_V</^ - VG't(y.). 
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Let ri^{t,x) = z^(t,x + c^(t)), then we have 

P-[v\to) + q ■ VGt,{v')] = g), 

P-[v\to) + q ■ VG,,{V')] = ^G,„(^i)(r,>, g), (6.25) 

Thus we obtain 

9Gt,{v'){^.q) - ^G,„(^i)(r,>,g) = <[<(to) + q ■ P-VG,„(^^)] 
= e'^^*»<^GoK(0))(r/^(0),^^(0)) - e^^*«</(0) • P_Vt;^(0) (6.26) 
+ <[rt{t,) - e^^*°77^(0)] + g • <P^VGtM)- 
The first term on the right of fl6.26p can be rewritten by using (16.161) 

<^G,{v-M,g-)=<&G,{v%){^% + nv',g') 

-<&G,{viM>, + nv\g') + &G,{v'){7^\<g-), 

where all functions are evaluated at t = 0, and the operator TZ is defined by 

n:=g^ ■ P>oVP^. (6.28) 

We say that a component in (I6.26P is negligible if its norm in E is much smaller 
than the right-hand side of (16.61) . So is the last term in (I6.26p . since A ^ 1 and 

||<P_VGi„(y.^)|U< p-i<Gi„((/.^)||« <m5^^ < (6.29) 

In order to estimate the other terms, we prepare rough bounds on the unknowns. 
Lemma [3.11 together with G G Qe.s implies that 



(6.30) 



ll^istr(0,l) < Wv'We, |I4o(0)IU < {l + GK + CS)\\v>o{to)\\E 

< {l + GK + G6 + Gi)yi,\\E. 
The estimates in (I3.12p together with G G Qe,s imply 

||<^^''||str(o,i) + ||<c'"||loo(o_i) < \\<w'"{to)\\n + \<c'"ito)\ 

< - n^ijE + Wy'llE + ^(w'>o + \b\My'\\E)) + \b\ (6.31) 



E) 



The equation for each {z^ ,g^) is given as follows. 

g = B', Z = F', 

X's := 26-\\\\wrj6'){w,z)^, I := {H{Q) - {V'Q.\w,))-\ 

B' := ;^^/ ^(^'^^^) + /((VQ,|zi) - g{V^Q,\w^))B + xsMz, Vw), ^3^32) 

, _( B'-VQ, + g-V^Q,-B 
■ \f'mzl-g■Vfmw^ + N' 

N' := x'sNc + XslWiQc + wi) - f'iQM^i ' 9 " VQJ + g ■ Vf'mw^] 
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where the superscript j and the dependence on {w, c) are omitted. Using the esti- 
mate ||ii'||str(-i,i) ^ ll'W^(^o) llw) we obtain in the same way as for fl3.12p . 

\x's\<\\4n/S, \B\<6', \B'\<6\\zU + \9\{\\wmnA6)\ 
\\F'h.n(-iA)<[T+{\\w{to)\\nAmMstr + \\9\\L^]+n^ 
In the case ||v5°||-h — ll'y^^'^ll'H ^ "we have = q, and so by the Strichartz estimate, 

\\^'\\str<\\z^ito)\\n+\q\\\w^{to)\\n- (6.34) 
In the other case Hv^'^Hh ^ ^, 

lk-''l|str+ II^^-'IIloo < \\z^{to)\\H + yito)\. (6.35) 
In both cases, f l6.4p implies 

Ik-'llstrCO,!) + ll5'''IU°°(0,l) 

< hh\\n + Wy^lu + + IqlMWv'h)] + \q\ (6.36) 

< \y>\\E+\q\My\\E). 



For the difference estimate, we consider the two cases separately. If Hv^'^H^ 
llyj^ll^ 3> S, then = q and 

<z = JV<z + (0, <f'{Qc>)zl - q ■ V/'(g)<< - q ■ <Vf'{Qc>)wl), (6.37) 

and so 

||<^1|str(0,l) < ll<^'(^0)|| + |&||k1|str+ |g|||<W^1|str+ |&r|g||k1|str, (6.38) 

where the last term comes from the last one of f l6.37p . Inserting f l6.3ip and fl6.36p . 
we obtain 

||<(2;^/)||strxL-(0,l) 

In the nonlinear case ||</3'^||^ ^ S, we have 

\<B'\ < ||<(w^c'')||^xR'^ll(^^^/^)lkxR'^ + lk1kll<(2;^/)|kxR'^, , . 

||<F'|Ui^ < T[||<^^||str + 5||<C^|l2°o|kl|L-] + lk1|str||<(;^^/)||strxL- 
+ ||<(W^ C") II strxL- II (z", /) llstrxL- , 

where the term with a power comes from the same term as in the linear case, i.e., 
(1 - X5)q ■ <Vf'{Qc>)wl. Hence 

||<(^",/)||strxL-(0,l) < ||<^"(to)||w+ [||<ti^1|str+ ||<C^|2d]||(^^/)||strxLo=, (6.41) 

which, together with fl6.3ip and (16.361) . leads to the same bound fl6.39p as in the 
linear case. 



34 K. NAKANISHI AND W. SCHLAG 

For the penultimate term in fl6.26p . we obtain from the equation of rj^ fl6.12p in 

the same way as above, 

\\<P-Kito)-e'^'°V-m\\E 

^ ll</l|L-||^^1lstr+ ||/||L-||m5t;1|L- (6.42) 

+ 5||<(^",/)||strxLr + ll<K,c")||strxLrll(-2''^?')llstrxLr' 

where we do not get the term with a-power, since the potential term is frozen in 
the 7] equation. Using f l6.30p -( l6.36p and Lemma 13.21 we can easily observe that the 
above is negligible because A ^ 1. Here again we have used that either Hv^'^H-h ^ ^ 
or = q, which will be tacitly utilized in the following, too. Hence the second term 
on the right of fl6.26p is also negligible by fl6.30p and (16.390 . 
For the remaining and leading term of fl6.26p . we have 

||gJOo^^G'o(^;")(r7^/)|U 

< e-'-'"A[\\vl, - r,,>4o|U + \<c^\MHoh)r'' ,^ 

><[\\vl, + nv'\\E + \g'\M\\vUE)] 

where t = 0. The penultimate term in (I6.43P is negligible thanks to £ <^ A, f l6.30p 
and fl6.39p . The last term in fl6.43p is dominated by 

\Kc>v' - Vik + II [P-. rlAAW + I<c1||7^^;l« + hUv^U 

< ||<^^||« + |<c^|||^°||^ + |<c^||^7°|||t;°||^ + I^^V^^'- 

Hence it is also negligible by using the estimates f l6.30p - fl6.39p and £ ^ A. 

It remains to deal with the leading term of f l6.43p , for which we need more precise 
estimates, employing the time decay of e"-*''. First we consider the linearized case 
||(/}°||-^ ~ llv^"'"ll'H ^ ^) using the equations 

^.(0) = e-^^*»G,„((^^), c°(t) = 0, c\t) = h, 

fj^ = JCr,^ + (0, g ■ Vf\Q)vi), g^{t) = q. 

The first component on the right of fl6.43p is estimated by 

II ^>o ^ T'lc^Ao We 

< ||e--^^*"P>o(v?° - nip')\\E + ||e-^^*«[P_,rb]<^ii5 + \\[e-''''^n]ipU\E (6.46) 
<e^^'^-ny^'\\E + C\b\y\\E, 
and the third component by 

\\ri% + nv^E < ||e"''^S>o(^o)IU + C|g|||v^°|U < ellV-IU + C|g|||v^°|U. (6.47) 

Since Hq ^ k and ^^dly^-'H^;) ^ H'/^-'ll^;, we see that the leading term of f l6.43p is 
smaller than 

^A - u^'We + |&|0.(||^li.)] [II^IU + kl0.(ll^'IU)] • (6-48) 

Hence the other terms, which have been shown to be negligible, are absorbed into 
the remaining half of f l6.6p . Thus the linearized case \\^p^\\ ^ 5 is done. 



(6.45) 



(6.51) 



(6.53) 
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It remains to consider the nonlinear case \\^p'^\\ ^ 6. Let 

e := w° - T^c>v^ = T*o<w'', (6.49) 

then we have 

^^>o - ^<c>^^>o = ^>o + [P-, r<ic'']v^, 

i = J£e + i?(^°,c°) ■ Ve + r;o[(0,</'(Q,.)W) +<^K,c^)]- ^^'^^^ 

The commutator term is neghgible thanks to f l6.30p . f l6.3ip and IIv^^Hk ^ ^- For the 
discrete spectral part of ^>o, we have 

\\Po4^ - e-'^*»e(to)]|U < 5||<(^^c^)l|strxLr(o,l), 
||e-^^%+(to)IU<e1|Po+e(to)IU, 
and for the continuous spectral part, 

dtic^lP.O/'^ = 5(^°,c°)[(/'(g)6|V6) - {mp.vo] 

+ {C^\P,T:om<fiQ,>)wl) + <Fiw'',cn]) (6.52) 

<5'llelll+llelU5||<«C^)||strxLr(0,l)- 

Thus we deduce 

c'^)||strxLi°=(0,l) 

<{1 + Ck + Cv^)[||¥^° - uv'We + l&l^^dlv^'IU)]- 
Similarly for the r]*^ component, fl6.12p implies 

h\h)\\E<\mE+c\q\\mE+cm\\E+\q\). 

||Po+[r/°(0) - e-^^*°r/0(to)]|U < 5|| (^°, ^?°) UstrxL^ , 
ai(£r/°|P7^°)/2 = B{w\c')[{VnQH\Vr^') - (£r^°|P,Vr^°)] 

+ (£P,r^°|(0,^?° ■ VnQH) + r;o[9(^,.)P(./;°,c°) ■ {z\g')]), 

and so, using ( 16.36^ we obtain 

||ry^o(0)IU <{1 + Ck + C^5 + Cl)\mE + C\q\{5 + I). (6.55) 

Also using fl6.30p - fl6.36p . we have 

\\nv'\\E<\g'\\\v'\\E<5m\E+\q\), 

05(|l4olU) <{1 + Ck + C5 + Cf.)U\W'\\E). 

|<c^(0)| < |6| + C\\<B{w\ &)\\Lr < n + - u^^We + l&l], 

|^?°| < |g| + C\\B^{w\ c')z' + B,{w\ c°)<7°|U^ < |g| + + |g|] 

Putting (16.531) . f l6.55p and the above estimates together, and using k + 
we see that the leading term in f l6.43p is bounded by 

- nv^E + 1610,(11^1,,)]"'' [II^IU + \q\mv'\\E)\ , (6.57) 

so the remaining half can absorb the negligible terms, concluding the proof in the 
nonlinear case Hv^'^IU ^ 5- □ 



(6.54) 



(6.56) 
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Mobile-differentiability of the fixed point G^, G Qi^s now follows from the closedness 
of Q^g for pointwise convergence. 

Lemma 6.4. Under the assumption of the above lemma, let Gn G ^"^^ be a sequence 
of maps such that Gn{'-p) G{(p) as n oo for all ip E H. Then G E G'^'s^- 

Proof. Since G„ E Qe^s, (16 ■4p implies that ^Gn{^p) for each (p E Ti is bounded in 
("H X M'^)*. Hence after extracting a subsequence, we have weak convergence of 
^Gnif) in ("H X R'^)* for (y9 in a dense countable subset A C Ti. To extend the 
convergence to all (p E H, take a sequence (pn converging to ip in "H. Then for any 
ip EH and q E M'', we have 

\\^Gk{'Pn){^,q) - ^Gk{^J{^,q)\\E 

< A||</.„-</.„r'H||^|U + |g|0,(||(/.„|U)] ^0, ^' 
as n, m — > oo, uniformly for all k E N. Hence we have the convergence 

lim ^Gki^p) = lim lim ^Gki'Pn) = lim lim ^Gk{(Pn) (Q.59) 

weakly in (T-L x W^)*. To see the mobile differentiability of G, we use the mean value 
theorem. For any (p,ilj E Ti, q E M.'^, k E N and £ G M small, there is 6' G [0, 1] such 
that 



Since Gk E QTt ■, we have 



G',.(<^(,)) - Gk{H>) = £^G',(/^)(^,g). (6.60) 

i,s . we ha 

||^G,.(/^)(^,g) - ^Gk{ip){ij,q)\\n 

<Ay''~vrn'mn + \q\M\\M]-^o 

as £ — )■ 0, uniformly for all k eN. Hence the limit G{ip) is mobile differentiable and 

^G{^){ij, q) = lim ^Gfc(^)(V', q), (6.62) 



fc— >CXD 

e,s ■ 

Therefore the fixed point belongs to Q"^'^ and so in particular C^'" in the Ti 



which implies G E QTt- ^ 



topology. Then it is easy to see that G^ is also C^'", so are A^c«,o and 

Appendix A. Table of Notation 





ordered pair and difference 


m 


(OA 5) 


minimum 






sum of two powers 


(16771) 


Br{0) 


a ball in the Banach space 




D,J,V 


basic operators 


(11.151). (12:71). Of 


f, d, p 


nonlinear ity, dimension and power 


dLH), dLZl) 


N{v),N{v),Nc{w), C{v) 


higher order nonlinearity around Q 


(12.201). (|2.2ip. Km 


A{v),Ac{w) 


transport terms 


(]2.26P, (13^ 


Xsiv) 


localizer around in 7^ 




M{v),Msiv) 


nonlinearity for v 


(|2.27|). ^ 
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B{w, c), F{w, c) 
U{t), U{t) 


nonlinearity for (c,w) 

localized flow in % and on graphs 


(I3.18|). Lemma Ea 


E{u),P{u) 
H, Tio, 'H± 
II • lb, K 

(•I-), '^(•'■) 

■ Str 


energy and momentum 
energy space and its subsets 
energy norm and the parameter 
bilinear forms on % 
Strichartz norm 


(|1.3|). (|2.22l) 

(OI). (|2.23l). (|2.24D 

(12.181). (12.171) 

(I23D, (USD, dMD 

(|3.10|) 


Q,Qc,Q{p, q), 

Q,Qip,q),y{Q),yo{Q) 
H{Q) = H^p{Q) 


the static solution, its transforms, 
vector forms, and the families 
kinetic energy matrix of Q 


(11.91). (11.121). (11.101) 
(|1.13|). (|1.14l). (|5.18[) 
(12. MP 


Pk,gk± 

Xk±,i^,iy,'y,v±,vd, . . . , 

P±,Po,Pd, P^ii P-y, ■ ■ ■ , 


linearized operators at Q 
eigenvalues and their bounds 
eigenfunctions of L+ , C 
spectral components of v 
the corresponding operators 


(\l.\b\],UM 
(|2.11D. 
(12.121). (|2.13l) 
(12.151). (12.161) 
(12.151). (|2.16p 


"10, ms, (j)s 
Tc, 


mobile distances and the cut-off 
translation operators 


(12.281). (12.381). (12.391) 
(I3.6P. (|3.28p 


\G\ 

Gif 

G 


sets of Lipschitz maps 
the graph in % 
the invariant map 
projection to "Hj, 


(I3.53P. (I3.69P. Ueiinition 16.21 
(|3.54p 

Theorem 13.61 
Lemma 15.11 


&G 

fje) 

Gt 

n 


mobile derivative 
translating variation in % 
complemented graphs 
spectral error of the (jr-translation 


Definition 16.11 
Definition 16.11 

dEED 

(I6.28P 
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